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Summary 

A set of conditions had not been formulated on the boundary of an elastic continuum since the time of 
Saint-Venant. This limitation prevented the formulation of a direct stress calculation method in elasticity 
for a continuum with a displacement boundary condition. The missed condition, referred to as the 
boundary compatibility condition, is now formulated in polar coordinates. The augmentation of the new 
condition completes the Beltrami-Michell formulation in polar coordinates. The completed formulation 
that includes equilibrium equations and a compatibility condition in the field as well as the traction and 
boundary compatibility condition is derived from the stationary condition of the variational functional of 
the integrated force method. The new method is illustrated by solving an example of a mixed boundary 
value problem for mechanical as well as thermal loads. 


Introduction 

The stress-strain law, the equilibrium equation (EE), and the compatibility condition (CC) are the 
three fundamental relations in elasticity. The material law was formulated in the mid-seventeenth century 
by Hooke (1635-1703). The equilibrium equation or the stress formulation is credited to Cauchy (1789— 
1857). Saint-Venant (1797-1886) developed the CC, or the strain formulation. It is a general belief that 
the fundamental elasticity relations were known for over a century. The thrust, therefore, was to develop 
approximate solution techniques because a closed-form solution cannot be generated for the vast majority 
of the solid mechanics problems. Such techniques included Airy’s method (ref. 1), Ritz’s method (ref. 2), 
the moment distribution technique (ref. 3), Kani’s method (ref. 4), the finite element technique (ref. 5), 
and others. 

It is surprising that the strain formulation was not known on the boundary of an elastic continuum, 
even though Cauchy’s stress formulation explicitly contained the boundary conditions also known as the 
traction conditions. Because of this deficiency, problems with displacement boundary conditions could 
not be solved using the direct stress calculation method, popularly referred to as the Beltrami-Michell 
formulation (BMF) (ref. 6). The strain formulation that was missed on an elastic boundary is referred to 
as the boundary compatibility condition (BCC). The BCC has been derived. Now the stress and strain 
formulations are parallel in form; each contain field equations and boundary conditions. Earlier, we 
derived the BCC for two-dimensional (ref. 7) and three-dimensional (ref. 8) problems in elasticity in 
Cartesian coordinates. The BMF was completed by adding the new BCC to the classical method. The 
completed Beltrami-Michell formulation (CBMF) can be used to solve displacement as well as mixed 
boundary value problems in elasticity. The CBMF stress formulation is as versatile as the Navier 
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displacement method, yet its equation structure is simpler. Solutions to plate and shell problems via the 
CBMF are discussed in references 9 to 11. 

A conservative elastician, believing the set of existing equations to be sufficient, may be reluctant to 
accept the new BCCs. However, it should be realized that some formulae and equations of the solid 
mechanics discipline were not completed in the first attempt, but were perfected eventually. For example, 
perfecting the flexure formulae required a century between Galileo, Bernoulli, and Coulomb. Saint- 
Venant completed the shear stress formula that was initiated by Navier. Cauchy formulated the stress 
equilibrium equation that was also developed by Navier in terms of displacement, but it contained only a 
single material constant instead of two. 

The formulation of the BCC in polar coordinates is the primary contribution of this report. The use of 
the new condition is also illustrated through the solution of a mixed boundary value problem for 
thermomechanical loads. The CBMF containing the BCC is obtained from the stationary condition of the 
variational functional n s of the integrated force method (IFM) (ref. 8). The variational calculation is 
performed in two distinct steps: 

(1) The terms of the functional tt v are transformed to obtain integrands, whose coefficients are either 
displacement variables, stress function, or reactions. 

(2) The stationary condition (ref. 12) of the functional 8n s with respect to displacement, stress 
function, and reaction yield all the expressions of the CBMF. 

The BCC is the coefficient of the variational stress function in the line integral term. Variational calculus 
in polar coordinates is more difficult than that in the Cartesian system (ref. 7) because the coefficients of 
the terms in the functional are functions of the r-coordinate. Also the Jacobian (J= r) has to be used. A 
nonvariational approach or carelessness can easily miss an expression because of the tensorial nature of 
stress and strain. The accuracy of CBMF derivation is essential because solution of elasticity problems in 
polar coordinates is very popular. Many existing elasticity solutions can be verified by back-substituting 
into the CBMF. To demonstrate the use of the new condition, two mixed boundary value problems are 
solved. The first example is for mechanical load, while the second is for thermal load. 

This report is organized as follows: First, a variational derivation is given for the CBMF. Green’s 
theorem is used for a quick validation of the new boundary condition. Then the CBMF is used to solve a 
problem with stress and displacement boundary conditions. This is followed by discussion and 
conclusions. Appendix A is a listing of symbols and acronyms found in this report. Appendix B presents 
the major steps of the variational derivation, which can be used by the reader to verify the BCC. 

Appendix C presents the solution of a structure using the IFM, which is the discrete analogue of the 
CBMF. 


Completed Beltrami-Michell Formulation in Polar Coordinates 

The CBMF in polar coordinates is obtained from the stationary condition of the variational functional 
(ref. 7) of the IFM. The functional n s has three terms (eq. (la)). The first term A (a, u ) represents the strain 
energy, expressed in terms of stress a and displacement u. The second term B(s, cp) is the complementary 
strain energy written in terms of the strain s and the stress function cp. The third term W is the potential of 
the work done. Basic steps of the derivation are given in appendix B. The functional is transformed into 
integrals with integrands whose coefficients are either displacement, stress function, or reaction variables. 
Symbolically it can be represented as follows: 


n s - A + B - W 


(la) 
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, = JJ* [ ( fi cld EE ) { w } + (field CC)(p]i/s + J [(boundary EE) {w}]jf + J[( boundary C-C- ) cp]c/ / 


D 


+ 


| (continuity condition) {reaction} d(. - 


= 0 


(lb) 


where the two displacement components j j are represented as {«}. Likewise, {reaction} represents the 
f reaction along r 1 

two reactions { >. Also, D is the plate domain; i x and are boundary segments where 

[reaction along 0J 

traction is prescribed and reaction is induced, respectively; and t is the line segment where stress is 
indeterminate. For a two-dimensional problem, £ = f’ 2 . The stationary condition of the functional in 
equation (lb) with respect to displacement, stress function, and reaction can be represented by the 
following symbolic expression: 


8n s =|| [(field EE) 8 {«} + ( field CC) Sep + | [(boundary EE)8{«}]e/f 
d q 

+| [(boundary CC) Scp]t/ f + j" (continuity condition) 8 {reaction} dt- 0 


(lc) 


The field EE and field CC are the coefficients of the variational displacement and stress function, 
respectively, in the surface integral terms of the functional (see also appendix B). Likewise the 
boundary EE and boundary CC are the coefficients of the variational displacement and stress function, 
respectively, in the line integral terms. The continuity conditions are the coefficients of the variational 
reactions. The field equations and boundary conditions of the CBMF recovered from the stationary 
condition of the variational functional in equation (lc) are as follows: 


Equilibrium equations 


Compatibility conditions 


Field: 


Field: 


do,. , 1 dx , K -tr e ) | b 
dr r 50 r 


= 0 


dx 1 0cr 0 

— + - 

dr r 00 


2x 

+ — + b e 


= 0 


Boundary: 


n r a r + n e x - P r 
n,.x + 77 0 a 0 = P Q 


(2a) 


(2b) 




1 0 2 a ; . 
r 2 00 2 


- u- 


0 2 a r 

dr 2 


+ - 


(1 + 2o) 0a ; . 
r dr 

(1 + 2u) 0G e 
dr z r 1 00 z r dr 

(l + u) d 2 x (l + u) dx 


8 2 <3a 


v 0 2 a 0 
r 2 00 2 


- + 


0r00 


dr 


= 0 


( 4 ) 


(3a) 

(3b) 


Boundary: 

0 / \ (l + u) f dx 

— (g 9 -og,.) — 

dr r \ 00 


+ g 0 - a r 


/, X0x , x da e 2(1 + 0 ) 

(l + o) (1 + O) -+ T 

V dr ' ’ 00 r 


n 0 =0 

( 5 ) 
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The kinematics boundary conditions are 


u-u =0 

(6a) 

c 

1 

Cl 

II 

o 

(6b) 


and 


r and 9 

are the polar coordinates 

a,., a 0 , and x 

are the stress components 

s,., Se, and y 

are the strain components 

u and o 

are the displacements 

b r and 6 0 

are body forces 

P r and P {) 

are tractions applied along boundary segment 1 1 

u and v 

are initial displacements along boundary segment 

9 

is the stress function 


The CBMF in equations (2) to (6) contains the following: 

(1) The stress formulation of Cauchy as coefficients of variational displacements 8 u and 8 u in the 
surface integral. It consists of two EEs in the field (eqs. (2a) and (2b)) and two on the boundary (eqs. (3a) 
and (3b)) that are popularly known as the traction conditions. The stress formulation has two distinct 
components: the field equations and the boundary conditions. 

(2) Strain formulation as the coefficient of variational stress function Sep. It is written in terms of 
stress for an isotropic material with Young’s modulus E and Poisson’s ratio u. It consists of the single 
field CC (eq. (4)), along with one new BCC (eq. (5)). Saint-Venant, unlike Cauchy, formulated only the 
field condition. Ele missed the boundary conditions that we have completed. Now, both the stress and 
strain formulations are consistent, containing the field equations (eqs. (2a), (2b), and (4)) as well as the 
boundary conditions (eqs. (3a), (3b), and (5)). In appendix B, the strain formulation is derived in terms of 
the strains. 

(3) Displacement boundary conditions. Two kinematics displacement boundary conditions (eqs. (6a) 
and (6b)) are obtained as coefficients of the variational reactions. A rigorous derivation of the continuity 
condition is more difficult than the stress and the strain formulations, which are straightforward. 

The three-component stress tensor (a r , cr 0 , and x) is indeterminate in the field and on the boundaiy 
because the state of equilibrium provides only two equations. To achieve determinacy of the stress state, 
we must add one CC in the field as well as one on the boundaiy. Saint-Venant has given us the field CC. 
We have formulated the BCC. For the derivation of elasticity equation, the variational technique is an 
elegant method because of the tensorial nature of stress and strain. A nonvariational approach may miss 
an equation. 

The CC should be imposed only when the domain is indeterminate, whether it is the field or the 
boundary. The CC has no relevance for a determinate domain or a determinate boundary. A plane stress 
problem is one degree indeterminate in the field because there are three stresses and two displacements. It 
has one field CC. A BCC should not be imposed on a free or a determinate boundary, where at least one 
stress component is zero. A clamped boundary is typically indeterminate; thus, one BCC is imposed. 

The solution of an elasticity problem using the CBMF has two distinct steps: 

(1) The stress state is calculated first using the equations (1) to (5) for an elastic continuum with stress 
and displacement boundaiy conditions. The displacement boundary conditions (eq. (6)) are not used, but 
the BCC is used. 
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(2) Displacements are back-calculated by integrating the strain field. The kinematics boundary 
conditions given by equation (6) are used to evaluate the constants of integration in the displacement 
function. 

The CBMF recognizes that displacement does not induce stress. The derivative of displacement, which 
becomes the strain that induces stress, is accounted for through the BCC. Displacement conditions are 
used to eliminate rigid body movement as explained in step (2). 

The BCC expressed in strain, stress, and displacement is as follows: 

Expressed in strain: 


e ; . 5 / n, 5y 

r rdr V 0 ' 2 r50 


n r + 


5s,. dy y 

rd 9 2 dr r 


h 9 =0 


(7a) 


Written in stress for an isotropic material with Poisson’s ratio u: 


d_ 

dr 


K 


- ua r ) 


(! + o)| 


5x 

1 - — W 

50 0 ' 


n„ + 


/i \ 5x / \ 5(7e 

v ’ dr v ’ 50 


2(1 + o) 

+ — -x 


h 0 = 0 (7b) 


In terms of displacements u and v: 


d 2 u r 5 2 u 
50 2 5r50 


5o 

50 


n r + 


d 2 u 2 

r l- u - r 

drdQ 


d 2 u 5o du 

— — - r r — 

dr 2 dr 50 


= 0 


(7c) 


The BCC expression contains either all three strains, or three stresses, or it contains derivatives of two 
displacement components. The BCC is not a continuity condition in displacement, stress, or strain; 
however, it is a function of the variables. As such, the BCC is expressed in the derivatives of stress, strain, 
and displacement, but it is not a component of rotation. The BCC is an independent condition. It forms a 
new elasticity expression that was missed since the time of Saint-Venant. The field CC is a second-order 
differential equation, while the boundary counterpart is a first-order equation. This characteristic is 
applicable to the stress formulation. The field EEs are first-order equations, while the boundary (or 
traction) conditions are algebraic equations. 


Annular Plate Subjected to Thermomechanical Load 


We will now illustrate the CBMF calculation strategy through the solution of a radially symmetrical 
annular plate with mixed boundaiy conditions for mechanical and thermal loads. Consider a plate made of 
an isotropic material with Young’s modulus E and Poisson’s ratio u. It has thickness h (considered unity) 
with outer and inner radii of a and b, respectively, as shown in figure 1(a). The mechanical load case 
consists of a uniform radial load of intensity p applied at the outer boundary r = a. The inner boundary is 
restrained: u = 0 at r = b. The CBMF for the mixed boundary value problem is generated from a special 
case of the variational functional. It is obtained using the condition of symmetry or by setting the shear 
stress x and transverse displacement u to zero (x = 0, u =0) as well as by neglecting variation with 


respect to the angle 0: 


5f . 

— = 0 . Also, a simpler stress function vj / is used. 
50 1 
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(a) Initial position under load p. (c) Displaced position 

(translation x and rotation 0). 


Figure 1. — Annular plate. 


b 

i 

s- 

ii 

s< 

b 

(8a) 

f d\\i 

CTn = r — — + — - U 
V dr r ) 

(8b) 


where U= 0 for this example. 


Solution for Mechanical Load Only 

The equations of CBMF for a symmetrical annular plate subjected to a uniform mechanical load of 
intensity p are listed in equations (9) to (1 1): 
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Equilibrium equations Compatibility conditions 


Field: 


Field: 


8a r | K-g 9 )_ o 
dr r 

(9a) 

d / n. (l + u L N. n 

(G 0 OG r )+ (G 0 G r ) = 0 

or r 

(10a) 

Boundary: 


Boundary: 


G,. - p at r = a 

(9b) 

g 0 - UG r = 0 at r = b 

(10b) 


where the displacement boundary condition is 

u = 0 at r=b (11) 

There are two stresses G,and Ge and one displacement u. In the field, there is one EE (eq. (9a)) and one 
CC (eq. (10a)). On the boundary there is one traction condition (eq. (9b)) and a single BCC (eq. (10b)). 
Also, at the inner boundary there is one kinematics condition (eq. (1 1)). Solution to equations (9a), (9b), 
(10a), and (10b) yield the stress response. The single displacement u is back-calculated by integrating the 
strain and evaluating the integration constant using the kinematics condition (eq. (11)). 

The EE and CC in the field are arranged to obtain the following two simpler uncoupled equations: 



+ g 0 ) = O 


(12a) 


da r | (o',- -cj 9 ) _ Q 

dr r 


(12b) 


Integration of the first equation (12a) yields the sum of the stresses a,- + Ge to be a constant. The second 
equation (12b) is uncoupled and solved. The two constants in the stress variables are determined from the 
traction condition (eq. (9b)) and the new BCC (eq. (10b)). The stress solution follows: 


_(l + u)r 2 +(l-u)h 2 ] f 2h 

(l + u)a 2 + (l - u)i 2 v r 2 y 


(13a) 


_(l + u)r 2 -(l-u)h 2 ] f 
(l + v)a 2 + (1 - v)b 2 yr 2 ^ 


(13b) 



(13c) 


(13d) 
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For a plate with E = 30x1 0 6 psi, o = 0.3, p = 1 psi, a = 20 in., and b = 10 in., the stresses at the outer 
boundary are a r = 1.0, while Ge = 0.763 psi and a, + Ge = 1.763 psi. At the inner boundary these are 
g, = 1.356 psi, Ge= 0.407 psi, and g,- + Ge = 1.763 psi. The sum of the stresses g, + Ge = 1.763 psi is 

62 

independent of the r coordinate of the plate. The BCC g 0 - ug,. = 0.62 — — has an inverse quadratic 

r 

variation with respect to the radius, with a minimum value of zero at the restrained boundary (r= 10 in.) 
and a maximum value of 0.46 psi at the outer free boundary (r = 20 in.). The BCC should not be imposed 
on the free boundary at r = a = 20 in. The stress state in the mixed boundary value problem is obtained 
without any use of the prescribed displacement boundary condition. The displacement function u is 
obtained following the standard elasticity solution strategy. Stress is changed into strain using Hooke’s 
law. It is integrated to obtain the displacement function that contained a constant c. 

u - 

£(l + u)a 2 + (l-u)6 2 

The constant c does not affect the stress state. For the problem, the constant is calculated to be zero 
(c = 0) from the homogeneous kinematics boundary condition u = 0 at r = b (eq. (11)): 

H 2 )b 2 -‘ 2 ) 

u — - 

E(l + v)a 2 +(l-u)b 2 

The values of displacements are u = 0 at the inner boundary, r = b, and u = 3.5x1 0“ 7 in. at the outer 
boundary, r = a. 

The CBMF produced the solution to the mixed boundary value problem in two steps: First the stress 
state was calculated using the field EE and CC, along with the traction condition as well as the BCC. 
Then the displacement function was back-calculated. Solution to the mixed boundary value problem 
could not have been obtained by the classical BMF stress formulation. Solution to the mixed boundary 
value problem is not available in standard textbooks in elasticity (refs. 6, 13, and 14). 


(14b) 


p + c 


(14a) 


Solution for Thermal Load Only 

The CBMF solution for the annular plate is obtained for a temperature distribution given as 

(X - T b ) 

T = T b + \ a ,\’ (r-b) 

(a -b) 


(15) 
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The temperature distribution is shown in figure 1(b). It has a linear variation with values T a and T h at 
r = a and r = 6, respectively. The coefficient of thermal expansion is a. The CBMF equations for the 
annular plate subjected to a thermal load are given below: 


Equilibrium equations 

Field: 

do,. | (g,. -c>9 
dr r 

Traction or boundary: 

a r = 0 at r-a (16b) 

where the displacement boundary condition is 


= 0 


(16a) 


Compatibility conditions 


Field: 

5 / x (1 + u ) / x 

— (a 9 - oa r ) + (a e - a r ) 

or r 


- -a E 


d T_ 
dr 


Boundary: 

a e - v a,. = -a ET at r - b 


u = 0 at r = b 


(17a) 


(17b) 


(18) 


Both the field and the boundary CCs (eqs. (10a) and (1 Ob)) for the mechanical load are modified for the 
temperature load to obtain equations (17a) and (17b). The field EE is not changed. The mechanical load is 
set to zero (p = 0) in the traction equation (16b). The EE and CC in the field are rearranged to obtain the 
following two simpler working equations: 


~-( a r +<J e)-~ aE ~r~ ( 19a ) 

dr dr 

du r | (a,. -cj e 
dr r 

The field equations are solved for the boundary conditions to obtain the response, consisting of a,., a () , and 
u: 


= 0 


(19b) 


a,. 


isa(a-r) 


3 r 2 ( a - b ) 


a 2 +b 1 + - b 2 j 


x| Tjt- 


+T b {~r 


a 1 +b 2 + o^n 2 - b 2 j 

+ rb 2 [n(l - u) -6(2 - o)] + ab 2 [n(l - o) - 6(2 - o)]| 

a" + b" + -b 2 j 

+ rb 2 [n( 2 + o) - 6(1 + o] + ab 2 [-6(1 + o) + a[2 + u)]j 


(20a) 
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E a 

3 r 2 {b-a) a 2 + 6 2 + o(n 2 - 6 2 ) 

x|V a j2r 3 a 2 +b 2 + v^a 2 -b 2 j - r 2 a 2 + 2b 3 + u(n 3 - 6 3 j + a 2 b 2 [a (l - u) - 6(2 - u)]| 

+T b |-2r 3 a 2 +b 2 + u^a 2 -6 2 j + r 2 a 3 - b 3 + 'iab" + v>[a 3 - b 3 ^ + a 2 b" []a(2 + u)-6(l + u)]jj 

(20b) 

a(l + u)(6 - r) 

3 r 2 (b-a) a 2 + b 2 + v^a 2 - b 2 jj 
x\T a |r 2 a 2 + b 2 + -b 2 j + ra 2 [a - 2b - u(a -6)] + a 2 6[a(l - u) -6(2 - u)]j (20c) 

+ T h |-7' 2 a 2 +b 2 + v(a 2 -6 2 j + ra 2 [2 a - 6 + 0(12-6)] + a 2 6[a(2 + u) - 6(1 + o)]jj 

The numerical values of the response parameters for T a = 100 °C, T b = 50 °C, and a = 12x10 6 /°C are 

(1) at r = a: a r = 0 ksi, g 6 = -17.5 ksi, and u = 0.012 in. 

(2) at r = 6: a, = 14.2 ksi, Ge = -13.7 ksi, and u = 0 in. 




The sum of the stresses a,- + Ge = 18508 - 1800r has a linear variation with respect to the r coordinate 
because of a similar distribution of temperature (see eq. (15)). The BCC 

1427.797 

g 0 - nc f + a ET — 6.478 + 0.780 r ksi is zero at the inner boundary. 

r 

The CBMF solved the thermal load problem for a mixed boundary value problem. Supeiposition of 
solutions for mechanical and thermal loads yields the result for thermomechanical combined load. 


Discussions 

This section examines the CBMF concept. Attributes of the CC are also given. The annular plate 
example is supplemented with an eight-bar discrete truss structure. The solution to the truss problem 
using the integrated force method (1FM), which is the discrete analogue of CBMF, is given in 
appendix C. 


Completed Beltrami-Michell Formulation 

Flooke’s law, which is common to all analysis methods, relates stress to strain through the material 
matrix [G]: 

{g}=[G]{s} (21) 
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Stress a must satisfy the state of equilibrium in the field as well as on the boundary of an elastic 
continuum. Likewise, strain has to comply with the condition of compatibility in the domain as well as on 
the boundary. The stress and strain formulations are sufficient for the determination of the stress state in 
an elastic continuum with stress and displacement boundaiy conditions. The equations that are required to 
calculate the stress state can be conceptualized in the following symbolic expression: 

Equilibrium equations 
Compatibility conditions 


{stress} ■ 


[ Mechanical load | 

I Initial deformation f 


( 22 ) 


The state of equilibrium and compatibility is sufficient for the determination of the stress state. 
Displacement is not required to calculate stress. An elastic body can undergo rigid body displacement and 
rotation that does not induce stress. Total displacement can be decomposed into an elastic component and 
a rigid body component: disp = disp e i as tic + disp n&A . The stress calculation in the CBMF accounts for the 
elastic component via the strain in the field and on the boundaiy. Recovery of the displacement from the 
stress state uses the kinematics or the rigid body displacement component. 

Calculating stress by combining the equilibrium and compatibility was envisioned by Michell, and it 
is described by Love in the following quotation (ref. 1): 

“It is possible by taking account of these relations [compatibility conditions] to obtain a complete 
system of equations [Equation 1 9] which must be satisfied by stress components, and thus the 
way is open for a direct determination of stress without the intermediate steps of forming and 
solving differential equations to determine the components of displacements.” 

The proposition of Beltrami and Michell can be realized now with the availability of the new BCC. A 
direct method is now available to calculate the stress state in a general elastic continuum with 
displacement as well as stress boundaiy conditions. The stress state is obtained without any recourse to 
displacement, which is back-calculated by integrating the strains. In the quotation, “intermediate steps” 
refers to Navier’s displacement method that contains higher order differential equations. For the annular 
plate example, the CBMF required the solution of two uncoupled differential equations. Navier’s method 
in contrast would have required the solution of a third-order differential equation. 

An 1FM of structural analysis has previously been formulated (ref. 8). 1FM is the discrete analogue of 
the CBMF in elasticity. In 1FM, forces are calculated from a set of equations [S]{.F} = { P* } that include 
the EE and the CC. Displacements are back-calculated. The 1FM solution to a truss problem is given in 
appendix C. 


Nature of Compatibility Condition 

The CC is a controller type of relation. Strains are controlled, f(s,., Se, y) = 0, in elasticity (or the 
strain formulation); likewise the deformations P are balanced, f(Pi, P 2 , . . . , P„) = 0, in a discrete 
structural system. The controller type of relation cannot be derived from an application of the standard 
concepts of mechanics, like “action equal to reaction” (leading to the EE), or the “cause effect relation” 
(that has given us Hooke’s law), or the “displacement continuity concept” (the “strain continuity” is 
conceptually incorrect). This is probably one important reason for the late development of these CCs. In 
elasticity, the field CC (or Saint-Venant's “strain formulation”) can be derived by simply eliminating 
the displacements from the strain displacement relations. However, the derivation of the BCC requires the 
use of variational calculus. For structures, a direct application of Saint-Venant's “strain formulation” 
would have been sufficient for the derivation of the CC (ref. 8). No calculus would have been required 
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because, like EE, the CC is also an algebraic equation. But such a procedure was not adopted, and the CC 
was not developed as a deformation balance concept. Variational calculus is the right tool to derive the 
BCC because of the tensorial nature of stress and strain. 


Nontriviality Property of Boundary Compatibility Condition 


The field CCs fee are satisfied automatically when expressed in continuous displacement functions u 
and o : fcc(u, v) = q(u, v ) - t(u, o ) = 0. However, the BCC, when expressed in terms of displacements, 
produces a nontrivial condition: 


' d 2 u r d 2 v 
v 50 2 drdQ 


dv 

56 


n r + 


d 2 u 

drdQ 


2 d 2 v dv du 
+ v — r — — - r r - 


dr 1 


dr 50 


«e = 0 


(23) 


The Navier displacement method should account for the BCC because this is not a trivial condition in 
displacement. The BCC should be enforced along the interelement boundaries in a finite element model. 
The role of BCC should be investigated further in the Navier displacement method. A two-span plate 
made of two different materials supported on an elastic foundation may be an ideal example for the 
investigation. 


Rotation and Compatibility Condition 


The BCC should not be confused with rotation. An elastic body under load moves from its initial 
position to occupy the final form by undergoing strain, an x translation, and a 0 rotation (see fig. 1(c)). 
Only strain (not translation or rotation) induces stress. Strain is zero when the body is rigid. Rotation and 
strain are independent of each other, even though both quantities are defined in terms of the derivatives of 
displacement. For example, the BCC, which is a function of the strains is defined in polar coordinates as 



8 ( \ 5y 

re a + — 

rdr V H ' 2r50 




5s,. 

r50 


_dy_ 
2 dr 


*e = o 


(24) 


The BCC enforces an equality constraint on the strain components; it imposes no restriction on either 
translation or rotation. The annular disk requires the BCC for analysis even while it is undergoing 
translation and rotation on a flat surface (with z = 0), as shown in figure 1(c). 
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Stability of Structure 


Consider the discrete truss shown in figure 2. Its analysis is given in appendix C. The truss has one 
field CC and one boundary CC as follows: 

Field CC: 


P 1+ P 2 -V2P 3 -V2P 4 + P 5 + P 6 = 0 (expressed in bar deformation, P) 


(25a) 


Oj + a 2 - 2o 3 -2o 4 +o 5 +o 6 =0 (expressed in bar stress, o) 


(25b) 


Boundary CC: 


P 2 + P 7 = 0 (in deformation) 

(25c) 

a 2 + a? = 0 (in stress) 

(25d) 


fy 


-►z 
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The field CC in equation (25a) restrains the six bar stresses. This is the discrete analogue of 
Saint-Venanf s strain formulation (see eq. (4)). On the boundary, two member stresses are related 
(eq. (25c)). The two stresses g 2 and g 2 cannot assume independent values in the lower boundary 
chord in figure 2(a). The situation is similar to the BCC of the annular plate. The stresses g,. 
and Ge cannot assume independent values along the inner boundaiy because of the BCC 
g 0 -u a,. = 0 => 0.4068 - 0.3 x 1.356 = 0, or g 6 = UG r . Compatibility conditions are required for the 
analysis of indeterminate structures, which are more stable than their determinate counterpart. Stability of 
a structure may concern the state of equilibrium. The CC does not degrade the stability of a structure. The 
original truss shown in figure 2(a) is displaced in two steps. The axial constraint at node 5 is released first. 
This process eliminates the BCC but the structure is still stable. It is one degree indeterminate and the 
support at node 5 can move along the x-coordinate direction. If the transverse restraint is released at that 
node, then the structure rotates. The truss undergoes an x translation and a 9 rotation as depicted in 
figure 2(b). It can be analyzed as a mechanism by accounting for the field CC (given by eq. (25a)). 


Number of Boundary Compatibility Conditions 

There are two displacement boundaiy conditions in a plane elasticity problem. The question is 
“should there be two BCCs?” The answer is “no.” The number of BCCs is equal to the number of field 
CCs, which is equal to the indeterminacy r, defined as the difference in the number of stress n and 
displacement m variables (r = n — m). The elasticity problem in polar coordinates has three stresses and 
two displacements. It is one degree indeterminate, and it has one field CC and one boundaiy CC. The 
number of displacement boundaiy conditions in the Navier displacement method is not equal to the 
number of BCCs. Consider the displaced position for the annular plate and the truss shown in figures 1(c) 
and 2(b), respectively. The plate undergoes translation and rotation, referred to as the kinematics 
conditions. The plate is restrained at the inner boundary, which is the elastic condition. The calculation of 
the displacement function in the Navier displacement method requires the simultaneous compliance of 
both types of conditions: kinematics as well as the elastic conditions. In the CBMF, the kinematics and 
the elastic conditions are enforced in two steps. First, the elastic condition, which essentially is the BCC, 
is used to calculate the stress response. In the annular plate example, the BCC is used to calculate stress 
in the CBMF. The kinematics condition u = 0 is then used to evaluate the integration constant in the 
displacement function. 


Concluding Remarks 

The boundary compatibility condition (BCC) for an elastic continuum has been derived in polar 
coordinates using a variational approach. The BCC in essence is a constraint that is imposed on the strain 
or the stress state. The new boundaiy condition completes the stress formulation in elasticity. The 
completed Beltrami-Michell stress formulation can be used to calculate the stress state in a general elastic 
continuum without any reference to the displacement in the field or on the boundary. The displacement is 
back-calculated from the stress state. The BCC when expressed in displacements yields a nontrivial 
condition. 
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Appendix A 
Symbols and Acronyms 


Symbols: 

A strain energy 

a, b plate outer and inner radii, respectively 

B complementary strain energy 

b r , be body forces 

c displacement function constant 

D plate domain 

E Young’s modulus 

f function 

{F} member force vector 

G r , Ge Green’s functions 

[G] material matrix 

h plate thickness 

J Jacobian 

£ domain boundary 

boundary segments 

n r , n e direction cosines 

{P*} load vector 

P r , Po prescribed loads or tractions 

p mechanical load 

R r , Re reactions 

r, z, 0 polar coordinates 

[5] IFM governing matrix 

T temperature distribution 
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U potential function 

u, v displacements 

u,v prescribed displacements 

V body force potential 

W potential of work done 

x translation 

a coefficient of thermal expansion 

P deformation 

s r , 80 , y plain strain components 
0 rotation 

4 function of displacement 

n s variational functional of the integrated force method 

a, bar stress 

g,., Ge, x plain stress components 

u Poisson’s ratio 

\| / simple stress function 

cp stress function 


Acronyms: 

BCC 

boundary compatibility condition 

BMF 

Beltrami-Michell formulation 

CBMF 

completed Beltrami-Michell formulation 

CC 

compatibility condition 

EE 

equilibrium equation 

IFM 

integrated force method 
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Appendix B 

Variational Formulation for the Completed Beltrami-Michell Formulation 

This appendix provides the variational derivation of the completed Beltrami-Michell formulation in 
polar coordinates that includes the new boundary compatibility condition (BCC). The equations are 
obtained from the stationary condition of the integrated force method (IFM) functional tc s , defined 
previously in equation (la) as 


7 i=A + B-W 


(Bl) 


where 


A = 


-# 


D 


du 0 x du S 0 x5o 0 Ga5u 0 


dr 


- + 


rdQ 


+ (JqU + 


dr 


■ + 


rdQ 


TO 

r 


[hrdrd 0 


(B2a) 


B = 


# 

D 


0 & r d 2 ip ® s ; .5(p 0 s 0 S 2 (p E0 yScp E3 y0 2 cp 


r 2 de 2 


+ 


rdr 


dr 1 


r 2 dQ rdrdQ 


hrdrdQ 


(B2b) 


W = h 


<6 


A 


uP r + o P Q 


j + ^> 

^2 


M 


uR r + o R e 


1161 

dt 2 + (U)( y b r ii + b e v) [rdrdQ) 


V 


D 


(B2c) 


The plate domain/) has boundary £, that is separated into segments and £ 2 ; i = l x + l 2 . Body 
forces are b r and be. Along the boundary segment £ ] , loads P, and P e are prescribed, and displacements 
wand o are free. The segment l 2 has prescribed displacements u and o that can induce reactions R x 
and R r . The derivation sets the uniform plate thickness to unity (h =1) without any consequence. 

The term A represents the strain energy, and it is expressed in stress and displacement, which are 
considered independent of each other. The strain energy term B is expressed in strain and stress function 
cp, which are also considered independent of each other. The potential of the work done is W. Body force 

dV dV 

potential V is defined as b r = and b {) = . The stress function cp is defined as 

dr rdQ 


= ap 

r rdr ,-ee 2 


(B3a) 


a s2(p V 
9 dr 2 


(B3b) 


T = — 


dr 


r dip N 
A'dQ y 


(B3c) 
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Each term of the functional is reduced to obtain new terms that contain two factors. The second factor 
can be displacement, a stress function, or reaction. The first factor is an expression in terms of stress, 
strain, and load. The stationary condition of the functional with respect to displacement, stress function, 
and reactions will yield the following expressions: 

(1) Field equilibrium equations (EEs) in stress. They are the coefficients of the variational 
displacements 8« and 8 o in the surface integral terms. 

(2) Boundary EEs, or traction conditions. They are the coefficients of 8 u and So in the line integral 
terms. 

(3) Field CC in strains. It is the coefficients of the variational stress function (Sep) in the surface 
integral term. 

(4) Boundary CC. It is the coefficient of (Sep) in the line integral term. 

(5) The displacement continuity condition. It is the coefficient of the variational reactions in the line 
integral term. 


Derivation of equations stated in items (1) to (4) listed above (see eq. (B4) below) is straightforward. The 
derivation of the continuity condition (item (5)) required back-calculation (see eqs. (BIO) to (B12) 
below). 

The first 1 1 terms of the functional reduced using techniques of calculus are given in equation (B4). 
The other five terms (12 through 16) are retained without any operation. 


0 
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(B4) 


All 1 1 terms are combined to obtain the following form of the functional: 
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The variation of the functional with respect to displacements 8 u and 8 u yields the field EEs 
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(B6a) 

(B6b) 
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Likewise, the field CC is obtained as the coefficient of the variation of the stress function Sep: 

S 2 &r de r | a 2 s e | 28e q dy d 2 y _ Q 

r 2 d 9 2 rdr dr 2 r Sr r 2 dQ rdrdQ 

Along the boundary segment i x , the variation of the displacements 8 u and 8 u yields the EEs or the 
traction conditions 


CT 77,. + T/7q =P r 


(B8a) 


T77 ; . + a 0 /7 0 = Pq 


(B8b) 


Along an indeterminate boundary, the BCC is obtained as the coefficient of the variation of the stress 
function Sep: 



d( rs e ) [ 

rdr 


2rdQ ) 


+ 


de r 
rd 9 


dy_ 

2dr 


«e = 0 


(B9) 


In summary, stress equilibrium is enforced in the field (eq. (B6)) and on the boundary (eq. (B8)). 
Likewise strain compliance is achieved in the field (eq. (B7)) and on the boundary (eq. (B9)). 


Displacement Continuity 

The displacement boundary conditions u-u - 0 and u - u = 0 are routinely used in analysis. Their 
derivations are shown through back-calculation. This strategy is followed to avoid artificiality in a direct 
derivation process. The expression (u - u ) 8 (P r = a r n r + x?7 0 ) + (u - o)S (P e = xn r + a e 77 0 ) = 0 yields 
the continuity conditions. Because IT and u are contained in terms 14 and 15 in equation (B2c), we have 
to prove the following formula along boundary segment i 2 '- 
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The variational form of equation (B 1 Oa) can be written as 
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Consider the reduction of the first of the two right-hand terms in equation B 1 Oa: 
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The terms b and c that contain higher derivatives of the stress function are reduced to obtain terms in the 
first derivative of the stress function: 
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The variation of the first right-hand term of equation (B 1 Oa) becomes 
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Likewise, the second right-hand term is reduced: 
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Verification of the formula given by equation (BlOa) is obtained by combining the two equations (Bile) 
and (B lid): 
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Verification of Boundary Compatibility Condition 

Green’s theorem is used for a quick verification of the BCC. The BCC is inserted in the line integral 
coefficient to recover the well known field CC in the surface integral term. The integral theorem in polar 
coordinates can be written as 


# (“ Jr( rG '- ) “ = § (°r n r + G e»e ) d( ’ 


(B13) 


where G, and Go are the coefficients of direction cosines n,- and ne in equation (B9), respectively: 
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(B14b) 
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The surface integral terms are generated as 
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(B15b) 


(B 16 ) 


The coefficient within the bracket is the field CC. The compatibility concept applies to the field as well as 
to the boundary. The nature of the compatibility expression changes in compliance with the domain and 
the boundary. The same interpretation is true for Cauchy’s field EEs. 
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Appendix C 

Solution to an Eight-Bar Truss 

The solution to the eight-bar truss shown in figure 2(a) is obtained using the integrated force method 
(IFM) (ref. 13), which is the discrete analogue of the completed Beltrami-Michell formulation (CBMF). 
The IFM, like the CBMF, generates the force solution by coupling the equilibrium equations (EEs) to the 
compatibility conditions. Displacements are back-calculated from the force solution. The truss is made of 
steel with Young’s modulus E = 30 000 ksi. Each of the eight bars has an area of 1 in 2 . Nodes 1 and 5 are 
fully restrained. It is subjected to a gravity load of magnitude P = -10 kip at the midspan location. The 
problem is to calculate the force and displacement response. 

The six EEs of the structure can be written in terms of bar forces F as 

00 0 ^ 1-100 0 
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10 0-^ 0 0 0 0 
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The two compatibility conditions (CCs) in bar deformations P can be written as 
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The CC is rewritten in member forces using the flexibility relation 



AE , 


for bar length L, area A, 


modulus E, and member force F: 
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Simultaneous solution of the six EEs and the two CCs yields the eight member forces { F } . The six 
displacements {X} are back-calculated: 
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